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Universal coarsening dynamics of a quenched ferromagnetic spin-1 condensate 
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We demonstrate that a quasi-two-dimensional spin-1 condensate quenched to a ferromagnetic phase under¬ 
goes universal coarsening in its late time dynamics. The quench can be implemented by a sudden change in 
the applied magnetic field and, depending on the final value, the ferromagnetic phase has easy-axis (Ising) or 
easy-plane (XY) symmetry, with different dynamical critical exponents. Our results for the easy-plane phase 
reveal a fractal domain structure and the crucial role of polar-core spin vortices in the coarsening dynamics. 


Ultra-cold atomic systems are well isolated from the en¬ 
vironment and present a pristine system for exploring non- 
dissipative manybody dynamics. An emerging area of explo¬ 
ration with these systems involves the dynamics induced by a 
quench across a phase transition to a symmetry-broken phase. 
Following the quench domains form, with each of these do¬ 
mains having made an independent choice for the symmetry¬ 
breaking order parameter. An aspect that has seen experi¬ 
mental investigation [1-4] involves quantifying the produc¬ 
tion of topological defects that emerge between domains im¬ 
mediately after the quench [5, 6]. Another aspect involves 
how these domains coarsen over time as the different broken- 
symmetry phases compete to select the equilibrium state. Of¬ 
ten at late times, when the domains are large compared to 
microscopic length scales, the coarsening dynamics is uni¬ 
versal: correlation functions of the order parameter collapse 
to a universal scaling function when the spatial coordinates 
are scaled by a characteristic length where t is the time 
after the quench [7]. The time dependence of this length 
scale L{t) ^ yields the dynamical critical exponent 
Most work in the classical theories of phase ordering kinet¬ 
ics has focussed on dissipative models relevant to temperature 
quenches. The late-time dynamics for systems undergoing 
conservative Hamiltonian evolution has developed as a new 
area of interest, particularly due to developments with ultra¬ 
cold atomic gases [8-13]. 

Spinor Bose-Einstein condensates exhibit both superfluid 
and magnetic order [14, 15] and present a rich phase diagram 
[16-18] for considering transitions between phases with dif¬ 
ferent symmetry properties. The simplest non-trivial case is 
a spin-1 condensate, which has been realised using ^^Rb and 
^^Na atoms with ferromagnetic and anti-ferromagnetic inter¬ 
actions, respectively. In general, an external magnetic fleld 
breaks the full spin symmetry of the Hamiltonian, reducing it 
to axial symmetry transverse to the fleld. The external fleld 
also leads to a quadratic Zeeman shift of the spin states that 
competes with the spin-dependent interaction to determine the 
equilibrium phase [16]. This system is ideally suited to study¬ 
ing phase transition dynamics because the Zeeman energy can 
be dynamically varied in experiments, allowing quenches be¬ 
tween phases [19, 20] (also see [21]), and because the sub¬ 
sequent magnetisation dynamics can be revealed with in situ 
imaging [22]. While the short-time dynamics following the 
quench is well-understood (e.g. see [23-27]), the subsequent 


domain coarsening has been the subject of considerable de¬ 
bate [23, 28, 29] and has been identifled as a signiflcant out¬ 
standing problem in the fleld [15]. 
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FIG. 1. (a) Phase diagram of a spin-1 BEC with < 0 and zero 
^-component of magnetisation indicating the transition between po¬ 
lar and ferromagnetic phases as the quadratic Zeeman energy q is 
varied. The quenches of q from qi to qf considered in this paper 
are indicated schematically, (b), (c) show the direction of the order 
parameter during late-time spin ordering for quenches into the easy- 
axis (b) and easy-plane (c) phases, with the color map indicated on 
the respective spin-spheres. Simulation parameters: no = 10^/^s, 
9n/\gs\ — 3. 

Here we study the non-dissipative dynamics of a quasi-two- 
dimensional ferromagnetic spin-1 condensate. We consider 
the quantum phase transition of this system from an unmagne¬ 
tised polar state to either an easy-axis or easy-plane ferromag¬ 
netic state, depending on the flnal value of the quadratic Zee- 
man energy (see Fig. 1). We demonstrate that both quenches 
behave universally in their late time coarsening dynamics. We 
And 1/z = 0.68 for the easy-axis phase, consistent with a bi- 
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nary fluid interpretation [30] and disagreeing with an earlier 
result of 1 / 2 ; ^ 1/3 [28]. A hydrodynamic analysis [31] also 
obtained a growth law, and showed that the growth re¬ 
duces to if the effects of superfluid flow are removed. A 
recent study of the coarsening dynamics of an immiscible bi¬ 
nary condensate revealed a growth law and verifled the 
scaling hypothesis by demonstrating correlation function col¬ 
lapse [13]. For the easy-plane case, we show that topological 
defects, namely polar-core spin vortices, play a crucial role 
in the dynamics, and And an exponent (accounting for loga¬ 
rithmic corrections to the coarsening) of ^ = 1.04, consis¬ 
tent with the z = 1 result for model E. An analysis of the 
structure factor for the quenches reveals an expected universal 
scaling with sl k~^ Porod tail for the easy-axis case, yet re¬ 
veals a non-integer tail in the easy-plane case. We verify that 
this arises from a fractal structure of the domains. Such frac¬ 
tal behaviour in the Porod tails has also been observed in the 
dynamical scaling of aggregates in dense colloidal solutions 

[32] . Our observation of coarsening in the easy-plane phase 
demonstrates that the system continually anneals towards an 
equilibrium state (c.f. [33]). 

The energy functional for a quasi-2D spin-1 condensate is 

[33] 




E= d^x 


+ ^\F? + 


( 1 ) 


where -0 = ( 7 / 1 , 7/05 The system has density inter¬ 

actions gnn‘^ (n = is the total areal density) and spin 
interactions gs\F\‘^ (F is the spin density with components 
Ffz = 'ipgfj.i’, where e {fxjyjz} are the spin-1 matri¬ 
ces). A magnetic held along 2 ; shifts the energies of the spin 
states. The linear Zeeman shift has been removed by trans¬ 
forming -0 into a frame rotating at the Larmor frequency. The 
quadratic Zeeman shift q can be tuned independently of the 
magnetic held using external microwave flelds (e.g. see [34]). 

For ferromagnetic interactions (gs < 0) the ground state of 
Eq. (1) can exist in three phases dependent on the relative val¬ 
ues of q and ngs [14]. The phase diagram for a system with 
zero magnetisation along 2 ; is shown in Fig. 1(a), and has been 
explored in experiments with ^^Rb [19, 29]. Here we consider 
the coarsening dynamics of this system quenched from a un¬ 
magnetized polar phase to a mdignQiizQd ferromagnetic phase 
by a sudden change in q. For 0 < q < 2|^s|no, where no is 
the initial condensate density, the magnetisation lies in the xy- 
plane (easy-plane) and we take f = {Fx, Fy)/no = F_\_/no 
as the order parameter. This phase breaks the continuous axial 
symmetry of the Hamiltonian, and the order parameter is not 
conserved. For g < 0 the magnetization lies along the 2 ; axis 
(easy-axis) and we take f = Fzjno sls the order parameter. 
This phase breaks the Z 2 symmetry of the Hamiltonian, but in 
this case the order parameter is conserved. 

To simulate the quench dynamics we numerically evolve 
the spin-1 Gross-Pitaevskii equations [14] with initial condi¬ 
tion of a polar condensate -0 = y^(0,1,0) that has vacuum 
noise added to Bogoliubov modes for q^ = 00 according to 
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FIG. 2. Results for easy-axis (left column) and easy-plane (right 
column) coarsening dynamics. (a,b) Growth in longitudinal and 
transverse magnetisation. Dashed lines indicate ground state lon¬ 
gitudinal (transverse) magnetisation for the easy-axis (easy-plane) 
cases. The excess energy per particle available for thermalisation 
is Qea = (|go - g/), Qep = |go(l - g//go)^, for the easy-axis 
and easy-plane quenches respectively, where go = 2|gs|no. For our 
parameters Qea ~ 4.5 x Qep, explaining greater thermal depletion 
in the easy-axis quench (c,d) Scaled order parameter correlation 
functions showing collapse, with unsealed data shown in the Insets. 
(e,f) Characteristic length scale L{t) extracted from G{r) (points) 
and best fits to results (lines). (g,h) Structure factor obtained from 
simulations scaled by the characteristic length scale (points). Best 
fits to high-Zc decay (lines). 

^ The initial noise in the condensate contributes an energy of approximately 
3 X 10“^Qea and 2 x 10 “^Qep in the EA and EP quenches, 
respectively. 
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the truncated Wigner prescription [33]. At t = 0 the quadratic 
Zeeman energy is suddenly reduced to its final value qf and 
dynamically unstable modes, seeded by the vacuum noise, 
cause magnetisation to grow exponentially [24, 33, 35] with a 
characteristic time scale tg = h/2\gs\no [Figs. 2(a), (b)] and 
characteristic domain size = ^/ \/2M|^s|no. The magne¬ 
tization growth saturates (at t ^ 10^ in our simulations) to¬ 
wards a value that is somewhat reduced from the ground state 
value due to a small thermal component that develops after 
the quench. This component arises from the thermalisation 
of the excess energy of the polar state over the ferromagnetic 
phase. As the magnetisation saturates we observe that the do¬ 
mains begin to coarsen and when they are large compared to 
^g [e.g. Figs. 1(b), (c)], the coarsening dynamics becomes uni¬ 
versal and independent of the microscopic details. 

In the universal regime we find, in agreement with the scal¬ 
ing hypothesis, that correlation functions of the order parame¬ 
ter have no explicit time-dependence when expressed in units 
of a characteristic length scale L{t) [7]. We examine the 
single-time correlation function 

G(r,i) = 2 J (fx{(l){x + r,t)-(l){x,t)), (2) 

for scalar or vector order parameter 0, where A is the area of 
the system and () denotes an ensemble average. We calcu¬ 
late this correlation function by averaging over an ensemble 
of 8 simulations that differ only by the random seeding. Our 
simulations are performed on grids with 1024 x 1024 (easy- 
axis) or 2048 x 2048 (easy-plane) points, using an adaptive 
step spectral method and we ensure that the simulations ac¬ 
curately conserve normalisation, energy and the ^-component 
of magnetisation. The long-time coarsening of the domains is 
revealed by the spreading of this correlation function, shown 
in the insets to Figs. 2(c) and (d). For the easy-axis phase we 
take L{t) to be the distance where G(r, t) first drops to zero, 
which is a measure of the average domain size. Because the 
easy-plane phase breaks a continuous symmetry, the notion of 
a single domain is not well defined. We therefore take L{t) 
to be the distance across which G(r, t) drops to 0.25G(0, t). 
Upon rescaling spatial coordinates by this length scale, the 
correlation function at different times collapse onto a single 
function G(r, t) ^ f as shown in Figs. 2(c) and 

(d), thus confirming the universal coarsening behaviour. 

We determine the dynamic critical exponent for the quench 
to the easy-axis phase by fitting to our results for L{t) 
[see Fig. 2(e)]. This yields 1/z = 0.68 in agreement with 
the growth law for a binary fiuid in the inertial hydro- 
dynamic regime [30] (also see [13]). We have found similar 
values for z (to within fitting errors) for simulations performed 
with Qf/\gs\no = { — 1-2, —1.8, —2.4}. The binary fiuid uni¬ 
versality class is also known as model H [36, 37]. 

For the easy-plane phase we determine the dynamic criti¬ 
cal exponent by fitting our results for L{t) to {t/ 

[see Fig. 2(e)]. This form has been used to describe coarsen¬ 
ing dynamics in the XY-model from an initial condition con¬ 
taining free vortices [38] (also see [39]), reflecting the slow 


approach to the asymptotic regime through the annihilation of 
vortex-antivortex pairs. In our simulations we observe that a 
large number of polar-core vortices [14, 40] emerge in the ini¬ 
tial unstable dynamics following the quench, which then de¬ 
cay away as 1/L{t)‘^ through vortex-antivortex annihilation. 
From this analysis we obtain 2 : = 1.04, which is consistent 
with z = 1 for model E in a two dimensional system [41] 
We have found similar values for 2 ; (to within fitting errors) 
for simulations performed with Qf /\gs\no = {0.2,1.2,1.8}. 
Model E describes a non-conserved planar ferromagnet dy¬ 
namically coupled to a second conserved field [36]. This 
fits our system well, where the second conserved field is Fz 
[23, 36, 42, 43]. Incorporating conservation of energy into 
model E gives model E', which also has 2 ; = 1 [44]. In the 
non-dissipative dynamics of a 2D (scalar) superfiuid a value 
of ^ = 1 was also observed [8]. 

It is also convenient to consider the order parameter struc¬ 
ture factor, obtained by Eourier transforming the correlation 
function 

S{k,t) = J d^rG{r,t)e^'^-^ = L‘^f{kL{t)), (3) 

where the scaling form follows from setting G{rA) = 
f {r/L(t)), with / being the Eourier transform of /. The 
structure factor is useful for examining the small r/L prop¬ 
erties of the order parameter, which can reveal the structure of 
domain walls and topological defects in the system [7]. Re¬ 
sults for the structure factor for the easy-axis and easy-plane 
quenches are shown in Eig. 2(g) and (h), respectively. 

Eor the easy-axis case we observe a “knee” in the structure 
factor at kL 1.3 followed by a “Porod tail” S{k) ^ k~^ 
for L > k~^ ^ ^g that indicates the presence of sharp do¬ 
main walls [7]. Eor small r/L, the probability of two points a 
distance r apart belonging to opposite domains is r/L, so that 
G(r, f) ^ 1 — 2r/L. This linear dependence on small r leads 
to the Porod law of decay for a d-dimensional system 

with domain walls. 

We also observe a Porod tail for the easy-plane case, but 
with a non-integer exponent, S ^ We interpret 

this non-integer Porod tail as arising from the domains hav¬ 
ing a fractal surface structure [45, 46]. Eor domains in a 
d-dimensional system having surface fractal dimension dg a 
j^-2d-\-ds emerges in the structure factor [47], reducing to 
the usual Porod law for the smooth surface case dg = d — 1. 
Thus our results in Pig. 2(h) suggest a surface fractal di¬ 
mension of dg ~ 1.5. To provide further evidence for this 
result, we determine a box-counting dimension for the do¬ 
main boundaries. We bin the easy-plane order parameter into 
discrete domains based on the spin direction and perform a 
box counting algorithm on the boundaries of these domains 


1 ^ 0.79 growth also fits the data in Fig. 2(f). Deviation between growth 

with and without a logarithmic correction only becomes apparent for 
t/ts > 10^, much longer than our simulations can investigate. 
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over an order of magnitude of box sizes. This yields a box 
counting dimension of D ^ 1.5 — 1.6. For comparison, we 
have also applied this analysis to the easy-axis domain bound¬ 
aries and extracted the box counting dimension ofD = 1.0. 
Possible physical implications of the fractal structure we ob¬ 
serve includes diffusion limited aggregation [48], or Schramm 
(stochastic)-Loewner evolution and the associated conformal 
invariance [49]. We note that the Porod tail in the easy-plane 
case is not accounted for by topological defects (vortices), 
which would result in a. k~^ tail [7, 50, 51]. 

A distinguishing feature of our system over more traditional 
models where coarsening dynamics has been observed is that 
our system has a firm microscopic foundation governed by 
conservative Hamiltonian evolution. It is therefore appeal¬ 
ing to explore the long-time microscopic details of our sys¬ 
tem dynamics. In particular, it is of interest to consider the 
role of vortices in the coarsening dynamics. We note that 
recently an exact analytic treatment of the vortex dynamics 
in the XY-model provided further insights into their role in 
coarsening dynamics [52]. Previous work has shown that the 
Landau damping rates for the spin-wave excitations is a slow 
and ineffective thermalization mechanism in the post-quench 
dynamics [33] but did not consider the role of spin vortices. 

Following the easy-plane quench we identify the decay of 
singly charged polar-core vortices. The state of a polar-core 
vortex is -0 ^ (sin^^e”*^, \/2 coS;5, sin^^e*^) where far 
from the vortex core cos (3 = ^/(l + q/2\gs\n)/2 [14]. The 
magnetisation lies in-plane with angle 0 that rotates by 27r/^ 
(k G Z) around the vortex centre, giving rise to a spin current 
but no mass current. At the vortex centre the particle den¬ 
sity concentrates in the 'ipo component (hence “polar-core”). 
While spin-1 condensates can support other vortices that com¬ 
bine mass and spin currents (e.g. Mermin-Ho vortices) [14], 
we only observe polar-core vortices of charge k. = ±1 (higher 
values of k, are unstable). The vortices are indicated in Fig. 3, 
revealing the decrease in vortex density as the coarsening pro¬ 
gresses. This occurs as (k, = 1) vortices and (k, = — 1) 
anti-vortices are drawn together and annihilate, leading to do¬ 
main annealing [see Fig. 3(c)]. The quantitative relationship 
between the vortex decay and the coarsening is revealed in 
Fig. 3(d). We also note that in the early stages of coarsening 
soliton like domain walls are observed as notches in the mag¬ 
nitude of the transverse magnetisation. These decay due to 
snake-like instability [25] that produces a (polar-core) vortex 
anti-vortex pair. A model for the interaction of spin vortices 
in a ferromagnetic condensate was proposed in Ref. [53], but 
here spin-waves appear to affect the dynamics. Thus a better 
understanding of the interaction of spin-waves with vortices 
is of interest, where it is possible that spin-waves provide an 
effective thermal field for the vortices during coarsening. 

In summary, our results for the Hamiltonian evolution of 
a spin-1 condensate quenched into a ferromagnetic phase re¬ 
veals a wealth of universal dynamics that can be controllably 
explored in experiments. Importantly, varying the final Zee- 
man energy changes the order parameter symmetry, hence 



x/^ x/^ 



FIG. 3. Polar-core vortices with = 1 (•) and k = —1 (0) in the 
(a) early stages and (b) later stages of the easy-plane coarsening dy¬ 
namics in a quadrant of the full simulation. Transverse magnetization 
indicated as in Fig. 1(c), but with saturation reduced to make vortices 
clear, (c) Evolution of the spatial region indicated with a dashed box 
in (b). The dashed boxes in (c) identify vortex-antivortex pairs that 
annihilate during the dynamics, (d) Total number of polar-core spin 
vortices (iVvort) as a function of time from a single quench simula¬ 
tion, demonstrating that the vortex density is proportional to L(t)“^. 
Other parameters as in Fig. 1(c). 


whether topological defects are supported, and also whether 
the order parameter is conserved during the dynamics: all cru¬ 
cial aspects of coarsening. The first steps towards coarsening 
dynamics in the regime we consider have been made in ex¬ 
periments [19, 35, 54], and will be aided by recent develop¬ 
ments of homogeneous trapping for cold gases [2, 3], which 
is advantageous for studying critical phenomena. An exciting 
aspect in experiments arises from the capability to observe 
vortices [1, 19, 55] and potentially track in situ domain dy¬ 
namics [56]. Dipole-dipole interactions will have a role in the 
dynamics of some spinor gases, although these can be elim¬ 
inated from the dynamics, e.g. using radio-frequency pulses 
[57]. Initial condition dynamics are known to be important in 
XY coarsening dynamics. The zero temperature quench we 
consider here only leads to the production of polar-core vor¬ 
tices, while Mermin-Ho vortices are expected in the case of 
a quench from a sufficiently high temperature initial condi- 
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tion and may change the coarsening dynamics [40] . Recently 
a suitable microscopic framework for simulating thermal dy¬ 
namics of spinor condensates has been developed [58]. 
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Zealand (contract number UOO1220). 


[1] C. N. Weiler, T. W. Neely, D. R. Scherer, A. S. Bradley, M. J. 
Davis, and B. R Anderson, Nature 455, 948 (2008). 

[2] N. Navon, A. L. Gaunt, R. P. Smith, and Z. Hadzibabic, Science 
347 , 167 (2015). 

[3] L. Chomaz, L. Gorman, T. Bienaime, R. Desbuquois, C. Weit- 
enberg, S. Nascimbene, J. Beugnon, and J. Dalibard, Nat Com- 
mun 6 (2015). 

[4] G. Lamporesi, S. Donadello, S. Serafini, F. Dalfovo, and 
G. Ferrari, Nat. Rhys. 9 , 656 (2013). 

[5] T. W. B. Kibble, Journal of Physics A: Mathematical and Gen¬ 
eral 9, 1387 (1976). 

[6] W. H. Zurek, Nature 317, 505 (1985). 

[7] A. Bray, Advances in Physics 43 , 357 (1994). 

[8] K. Damle, S. N. Majumdar, and S. Sachdev, Phys. Rev. A 54 , 
5037 (1996). 

[9] B. Zheng, Phys. Rev. E 61 , 153 (2000). 

[10] K. Koo, W. Back, B. Kim, and S. J. Lee, J. Korean Phys. Soc. 
49 , 1977 (2006). 

[11] A. Asad and B. Zheng, J. Phys. A 40, 9957 (2007). 

[12] H. Takeuchi, K. Kasamatsu, M. Tsubota, and M. Nitta, Phys. 
Rev. Lett. 109 , 245301 (2012). 

[13] J. Hofmann, S. S. Natu, and S. Das Sarma, Phys. Rev. Lett. 
113 , 095702 (2014). 

[14] Y. Kawaguchi and M. Ueda, Physics Reports 520, 253 (2012). 

[15] D. M. Stamper-Kurn and M. Ueda, Rev. Mod. Phys. 85 , 1191 
(2013). 

[16] J. Stenger, S. Inouye, D. M. Stamper-Kum, H.-J. Miesner, A. P. 
Chikkatur, and W. Ketterle, Nature 396, 345 (1998). 

[17] T.-L. Ho, Phys. Rev. Lett. 81, 742 (1998). 

[18] T. Ohmi and K. Machida, J. Phys. Soc. Jpn 67, 1822 (1998). 

[19] L. E. Sadler, J. M. Higbie, S. R. Leslie, M. Vengalattore, and 
D. M. Stamper-Kurn, Nature 443, 312 (2006). 

[20] E. M. Bookjans, A. Vinit, and C. Raman, Phys. Rev. Lett. 107 , 
195306 (2011). 

[21] S. De, D. L. Campbell, R. M. Price, A. Putra, B. M. Anderson, 
and 1. B. Spielman, Phys. Rev. A 89, 033631 (2014). 

[22] J. M. Higbie, L. E. Sadler, S. Inouye, A. P Chikkatur, S. R. 
Leslie, K. L. Moore, V. Savalli, and D. M. Stamper-Kurn, Phys. 
Rev. Lett. 95 , 050401 (2005). 

[23] A. Lamacraft, Phys. Rev. Lett. 98, 160404 (2007). 

[24] H. Saito, Y. Kawaguchi, and M. Ueda, Phys. Rev. A 76,043613 
(2007). 

[25] H. Saito, Y. Kawaguchi, and M. Ueda, Phys. Rev. A 75,013621 


(2007). 

[26] M. Uhlmann, R. Schiitzhold, and U. R. Fischer, Phys. Rev. 
Lett. 99, 120407 (2007). 

[27] B. Damski and W. H. Zurek, Phys. Rev. Lett. 99, 130402 
(2007). 

[28] S. Mukerjee, C. Xu, and J. E. Moore, Phys. Rev. B 76, 104519 
(2007). 

[29] J. Guzman, G.-B. Jo, A. N. Wenz, K. W. Murch, C. K. Thomas, 
and D. M. Stamper-Kurn, Phys. Rev. A 84, 063625 (2011). 

[30] H. Furukawa, Phys. Rev. A 31, 1103 (1985). 

[31] K. Kudo and Y. Kawaguchi, Phys. Rev. A 88, 013630 (2013). 

[32] T. Sintes, R. Toral, and A. Chakrabarti, Phys. Rev. E 50 , R3330 
(1994). 

[33] R. Barnett, A. Polkovnikov, and M. Vengalattore, Phys. Rev. A 
84 , 023606 (2011). 

[34] E. Gerbier, A. Widera, S. Eolling, O. Mandel, and I. Bloch, 
Phys. Rev. A 73 , 041602 (2006). 

[35] S. R. Leslie, J. Guzman, M. Vengalattore, J. D. Sau, M. L. Co¬ 
hen, and D. M. Stamper-Kum, Phys. Rev. A 79,043631 (2009). 

[36] P. C. Hohenberg and B. 1. Halperin, Rev. Mod. Phys. 49, 435 
(1977). 

[37] A. Bray, Phil. Trans. R. Soc. A 361, 781 (2003). 

[38] A. J. Bray, A. J. Briant, and D. K. Jervis, Phys. Rev. Lett. 84 , 
1503 (2000). 

[39] B. Yurke, A. N. Pargellis, T. Kovacs, and D. A. Huse, Phys. 
Rev. E 47 , 1525 (1993). 

[40] K. Kudo and Y. Kawaguchi, Phys. Rev. A 91, 053609 (2015). 

[41] B. 1. Halperin, P. C. Hohenberg, and E. D. Siggia, Phys. Rev. 
Lett. 32, 1289 (1974). 

[42] K. Nam, B. Kim, and S. J. Lee, J. Stat. Mech 2011, P03013 

( 2011 ). 

[43] U. Tauber, Critical dynamics (Cambridge University Press, 
2014). 

[44] B. 1. Halperin, P. C. Hohenberg, and E. D. Siggia, Phys. Rev. 
B 13, 1299 (1976). 

[45] H. D. Bale and P. W. Schmidt, Phys. Rev. Lett. 53, 596 (1984). 

[46] D. W. Schaefer and K. D. Keefer, Phys. Rev. Lett. 56, 2199 
(1986). 

[47] C. Oh and C. Sorensen, J. Nano. Res. 1, 369 (1999). 

[48] T. C. Halsey, Phys. Today 53 , 36 (2000). 

[49] J. Cardy, Ann. Phys. 318 , 81 (2005). 

[50] A. J. Bray and S. Puri, Phys. Rev. Lett. 67 , 2670 (1991). 

[51] F. Rojas and A. D. Rutenberg, Phys. Rev. E 60, 212 (1999). 

[52] A. Forrester, H.-C. Chu, and G. A. Williams, Phys. Rev. Lett. 
110 , 165303 (2013). 

[53] A. M. Turner, Phys. Rev. Lett. 103, 080603 (2009). 

[54] M. Vengalattore, J. Guzman, S. R. Leslie, F. Serwane, and 
D. M. Stamper-Kurn, Phys. Rev. A 81, 053612 (2010). 

[55] S. W. Seo, S. Kang, W. J. Kwon, and Y.-i. Shin, Phys. Rev. 
Lett. 115 , 015301 (2015). 

[56] K. E. Wilson, Z. L. Newman, J. D. Lowney, and B. P. Anderson, 
Phys. Rev. A 91, 023621 (2015). 

[57] M. Vengalattore, S. R. Leslie, J. Guzman, and D. M. Stamper- 
Kurn, Phys. Rev. Lett. 100 , 170403 (2008). 

[58] A. S. Bradley and P. B. Blakie, Phys. Rev. A 90, 023631 (2014). 



